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Abstract 

The induced charge density, Pi n d{r), generated in graphene by the potential well of the finite 
radius R is considered. The result for Pindir) is derived for large distances r S> R- We also 
obtained the induced charges outside of the radius r>fi and inside of this radius for subcritical 
and supercritical regimes. The consideration is based on the convenient representation of the 
induced charge density via the Green's function of electron in the field. 
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I. INTRODUCTION 



As known, the induced charge density, Pi n d(r), in the external electric field appears due 
to vacuum polarization. In the field of heavy nucleus, this important effect of Quantum 
Electrodynamics (QED) was investigated in detail in many papers, see, e.g., Refs.-~-. New 
possibilities to study vacuum polarization in QED at large coupling constant have appeared 
after recent successful fabrication of a monolayer graphite (graphene), see Ref.- and recent 
Review^. The single electron dynamics in graphene is described by a massless two-component 
Dirac equation^— so that graphene represents a two-dimensional (2D) version of massless 
QED. On the one hand, this version is essentially simpler than conventional QED because 
effects of retardation are absent due to instant Coulomb interaction between electrons. On 
the other hand, the "fine structure constant" a = e 2 /hvp is of order of unity since the 
Fermi velocity vp ps 10 6 m/s ~ c/300 (c is the velocity of light), and therefore we have a 
strong-coupling version of QED. Below we set h — c — 1. 

Screening of charged impurity in graphene can also be treated in terms of vacuum 
polarization^— . Investigation of impurity screening is important for understanding of the 
dependence of transport properties on the impurity concentration. There are two regimes for 
the Coulomb impurity in the gapless graphene, subcritical and supercritical. In the subcrit- 
ical regime, it is shown in the leading order in a and exactly in the Coulomb potential that 
the induced charge is localized at the impurity position, see^r— >2&2i. In the supercritical 
regime, vacuum polarization in the Coulomb field has been recently considered in Refs.— 
In this case, the induced charge density is not localized at the impurity position due to the 
effect similar to that of e + e~ pair production in 3D QED in the electric field of supercritical 
heavy nuclei. In the present paper, we answer to the question whether the phenomenon of 
the induced charge localization also exist in the potential well of finite size R and depth Uq. 
Namely, we calculate the asymptotics of Pi n d{v) in the field of an azimuthally symmetric 
potential well at large distances r 3> R. We apply the method suggested in Ref.— for calcu- 
lation of the finite nuclear size effect on the induced charge density at large distances in a 
strong Coulomb field in 3D QED. We show that there are also subcritical and supercritical 
regimes in this problem. However, the induced charge is not localized at r < R in the 
subcritical regime and has power "tail" in contrast to the case of the Coulomb field. In the 
vicinity of transition from the subcritical regime to the supercritical one, small variation of 
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the potential parameters drastically changes the induced charge density. We demonstrate 
that this fact is not related to the smoothness of the potential well. The attempt to calculate 
the induced charge distribution in the potential well in graphene was previously performed 
in Ref.—. The authors of this paper used the method which akin to that used at calculation 
of conventional Friedel oscillations. However, our results for the induced charge density 
differ from that obtained in Ref.— mainly due to the mistake performed in Ref.— at the 
calculation of the phase shift. 

The paper is organized as follows. In Section [III we derive the general expression for 
the induced charge density convenient for calculation of the asymptotics at large distances. 
In Section II I II we consider the Green's function of electron in an azimuthally symmetric 
potential and use this function in calculations of Pi n d(f) in Section HVl Critical values of g 
are discussed in Section [V] calculating the scattering phase shifts of electron wave function 
in the field of the potential well. The induced charges outside of the radius r ^> R and 
inside of this radius for subcritical and supercritical regimes are considered in Section IVIl 
Finally, in Section IVIII the main conclusions of the paper are presented. 

II. GENERAL DISCUSSION 

In graphene, the induced charge density in the potential U(r) have the form 



where N = 4 reflects the spin and valley degeneracies, and the Green's function G(r,r'\e) 
satisfies the equation 



Here a = (eri,<72), and <7j are the Pauli matrices; p = (p x ,Py) is the momentum operator, 
r = (x, y), and I = diag{l, 1}. The matrixes er act on the pseudo-spin variables and the spin 
degrees of freedom are taken into account in a factor N. According to the Feynman rules, 
the contour of integration over e goes below the real axis in the left half-plane and above 
the real axis in the right half-plane of the complex e plane. Using the analytical properties 
of the Green's function, we deform the contour of integration with respect to e so that it 
coincides finally with the imaginary axis. Then we follow Ref.— and write the equation for 




(1) 



[e — U (r) — vf<t ■ p] G(r, r'\e 



) = 6(r- r')I. 



(2) 
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the Green's function in the form 



G{r,r'\ie) = G (0 \r,r'\ie) + J d ri dr 2 G {0 \r , r\\ie) [U(ri)5(ri — r 2 ) 

+ U( ri )G(r u r 2 \ie)U(r 2 )} G^(r 2 , r'\ie) , (3) 



where G^°\r, r'\ie) is the solution of Eq.(j2J) at zero external field. 
It is convenient to represent Pindiy) as a sum 



PiUr) = P&(r) + p&(r) , (4) 

where p\nd( r ) * s the hnear in U(r) contribution and p^ d {r) is the contribution of high order 
in U(r) terms. It follows from Eqs.([T]) and ([3]) that 



Pind( r ) 



eN J°° ^-J dr{Ir^G^\r 1 r 1 \ie)U{r l )G^{r ll r\ie)Y (5) 



Formulas (JH]) and ([6]) are very convenient for calculation of the induced charge density at 
large distances. 



III. GREEN'S FUNCTION IN AN AZIMUTHALLY SYMMETRIC POTENTIAL 

The free Green's function G^°\r, r'\ie) is given by 

(rr . t\ 1 

(7) 



G(°>(r,r'|ze) = 



K Q {\e\i)-^{e)^^K x {\e\i) 



2vr 

where £ = r — r', and K Qt i(x) are the modified Bessel functions of the third kind. Let 
us represent the electron Green's function G(r,r'\e) in an azimuthally symmetric potential 
U (r) in the form 

1 ^ ■ ^ I A m (r,r'\e) -ie'^' B m (r, r'\e)\ 
G(r,r'\e) = — V e™**-*^ , (8) 

ZkJT'oo \ie^C m (r,r'\e) e^-^AJr, rV 

and use the relation 



S( r - r') = 5{r y e im ^ . (9) 

27rvrr / 



OG 



m=— oo 
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Then, from Eq.(j2]) we obtain the equations 



u(r))Am -?^.^±lc m i{r - r,) 



dr r 



(e-U(r))C m + ^--A m = 0, (10) 
Or r 

and the relations D m = A_ m _i and B m = — C_ m _i. Therefore, to find the Green's function 
in a azimuthally symmetric potential, it is sufficiently to solve equations (fTDl) . 



IV. AN INDUCED CHARGE DENSITY AT LARGE DISTANCES 

To calculate the asymptotics of the function p\^ d at distances r>i?, where R is a typical 
size of the potential, we can put r 1 = in the arguments of the free Green's functions in 
Eq-©. After that we take the integral over e and obtain: 

P&(0 = ^ f dr'r'U(r'). (11) 



16 r 3 

One can see that the induced charge density in the leading order in the external field goes 
to zero at large distances as 1/r 3 . 

Let us consider the function p^ d {r) at r R. We substitute Eqs.([H]) and ((7]) to Eq.©, 
put ri = and 7*2 = in the arguments of the free Green's function, and take the integral 
over angels of the vectors t*i and r 2 . Then we obtain 
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PiOd(r) = -^ J dee 2 [K 2 (\e\r) - K 2 (\e\r)] J J dridr 2 r 1 r 2 U(r 1 )U(r 2 )A (r 1 ,r 2 \ie) . (12) 
-00 

Here Ao(r\,r 2 \ie) is the solution of Eq.f llOp at m — 0. Note that Eq.f lT2|) includes the 

contributions of the terms with m = and m = — 1 in the Green's function ([8]), since 

D_i = A . It is convenient to introduce the functions 

a(r, e)= / dr'r'U(r')A (r, r'\ie) , c(r, e) = / dr'r'U(r')Co(r, r'\ie) , (13) 
Jo Jo 

Let us multiply both sides of the equations fllOp by r'U (r'), and then take the integral over r' 
from zero to infinity As a result we obtain the following equations for the functions a(r, ie) 
and c(r, ie): 

c tt( \\ t \ (5c ( r ' e ) c ( r ' e ) TT( \ 
(ie - U (r))a(r, e) = U{r), 

or r 

(i e _ U(r))c(r, e) + = . (14) 
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The boundary conditions for these equations are a(0, e), c(0, e) < oo, and lim a(r, e) = 
lim c(r, e) = 0. In terms of the function a(r, ie), Eq. ffl2l) has the form 

T— >OQ 

AT /"OO /"OO 

Pfl(r) = d" 2 [Kli\e\r) - tf 2 (|e|r)] ^ *V(/(r')a(r', e) . (15) 

Then we pass in this equation from the variable e to the variable E = re and replace 
a(r',E/r) on a(r',0) at r ^> R. We can do that because the integral over E converges at 
E ~ 1 due to the properties of the If- functions. After this replacement we take the integral 
over E and arrive at the following expression for the asymptotics of p ind (r): 

AT POO 

P { 2 d (r) = — / dr'r'U(r')a(r\0) . (16) 



16 r 3 j 

Thus, the function p^ d {r) has the same behavior at large distances as p^Jj)- 

Let us consider a simple example of the potential, U(r) = —Uq6(R — r), where 9(x) is 
the step function, R is the radius of the potential well. The solution a(r, 0) of Eq. (fl4|) is 



a(r,0)={ J ^ ' ' (17) 




where J n (a:) is the Bessel function. Using this solution, we find the sum of the contributions 
Eq-ffTTD and Eq. flTB]) . 

where g = UqR is the effective dimentionless coupling constant. The induced charge density 
( fl8l) is the odd function of the parameter g, which corresponds to the Furry theorem in 
QED. The formula ( Tl8l) contains singularities at the critical values of g = g c satisfying the 
equation Jo(g c ) = 0. In our case, the first three values are g c ~ 2.41, 5.52, 8.65. Existence of 
such singularities is not related to strong variation of our potential around the point r = R. 
We found numerically the first three critical values of g = UqR for the smooth potentials 
U^r) = -U e- r l R and U 2 (r) = -U e- r ^ R2 . In the first case, g c « 2.87, 5.9, 9.0. In the 
second case, g c ~ 2.7, 5.7, 8.0. We see that the corresponding numerical values of g c are 
close to each other. 

Actually, singularities in Eq. (I18j) have appeared as a result of substitution a(r',E/r) — > 
a(r', 0) in Eq. (fl6|) . which is not valid in the vicinity of g = g c since Eq.( fiT|) has no sense at 
g = g c . In the vicinity of g = g c , it is necessary to perform calculation of the integrals in 

6 



Eq.( !T5|) more accurately. For the step-like potential U{r) = —Uq6{R — r), the solution of 
Eq. fTHj) at e 7^ has the form 

, lj (( Uo + ie)r)-U /(U + ie), r<R, 
a(r,e) = { 7 (19) 
PK (\e\r), r>R, 

, ±Ji((£/" + z'e)r) , r < R, 
c( r ,e) = <( 7 m ; ; (20) 

-i/3sign(e)i^i(|e|r) , r > R. 



Taking into account continuity of the functions a(r, e) and c(r, e) at r = i?, we obtain 
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J ((t/ + ie)R) - iBign(e) Ji((i7o + 



(21) 



^i(l^). 

Then we substitute Eqs.f ll9p and (}2~Tj) to Eq. (fT5|) and take the integral over r'. As above, 
the main contribution to the integral over e at r 3> R is given by the region e < 1/r, so that 
we can use the relations eR <C 1 and e/[/o <C 1- Finally we find the expression for the sum 
of p\n d (r) and p\nd( r ) a ^ l ar g e distances, 

_ eNMg)Mg)R [°° , jg(er) - if 2 (er) 

This expression is valid at arbitrary value of the coupling constant g = UqR. If |Jo(#)| ^ 
(R/r)\n(r / R), then it is possible to neglect the second term in the denominator of the 
integrand, and we return to the expression ( fl8l) . If g is close to some g c so that | Jo(g)\ 
(R/r) \n(r/R) <C 1, we obtain 

eiVsign (g - g c ) 

pn r ) = ~o — r 23 

27rr / In \g — g c \ 

In this case the induced charge density diminishes as 1/r 2 and has opposite sign for g < g c 
and g > g c . In terms of distances, the asymptotics (|23|) is valid at 

1 < r/R<&-hx\g-g c \/\g-g c \ . 

At r/R^$> — In \g — g c \/\g — g c \ and \g — g c \ <C 1, we have (see Eq. (TT8Tl ) 

eNR 

»v = m^:y < 24 > 

In order to illustrate the transition from the asymptotics ( 1231) to the asymptotics (l24p . we 
consider the ratio Pi n d( r ) I Piif) at r 3> R and \g — g c \ <C 1. In this case this ration depends 
only on the variable rj = —r\g — g c \/(R\a \g — gc\)- The dependence of Pi n d(r) / p2(r) on 77 is 
shown in FigJTJ We see that Pi„d(r) ~ p-iir) already at g ~ 2. 




FIG. 1: Ratio Pind(f) / P2(f) at r 3> R and \g— g c \ <C 1 as a function of 77 
The asymptotics p2(r) is given by Eq. ([M|) and pi n d by Eq. ([22j) . 



-r\g-g c \/(Rln\g-g c \). 



V. CRITICAL VALUES OF g AND SCATTERING PROBLEM 

It is possible to explain critical values of g using the approach based on the scattering 
problem, as it is usually performed at the consideration of Friedel oscillations, see Ref.— . 
Writing the wave function of electron as 



ij)(r) 



u m (r)e im * 
%d m (r)e^ m+l ^ 



(25) 



we obtain equations for the functions u m (r) and d m (r), cf. Eq.(TjJJ 

dd m m + 1 



U(r))u 
U{r))d, 



m n 

or 
du ri 



r 



dm ! 



m 

-7; u m = 

or r 



(26) 



The solution of this equations in the the step-like potential has the form (common nor- 
malization factor is omitted): 



u m (r) 



J m (\U + e\r), r<R, 
^mJm{W) + y m N m (\e\r) , r > R, 



(27) 



d m (r) 



sign([/ + e) Jm+iG^o + e\r) 



r < R 



(28) 



sign(e) [pL m J m+1 (\e\r) + u m N m+1 (\e\r)] , r > R. 
Here N m (x) are the Bessel functions of the second kind. From continuity of the functions 



S 



u m {r) and d m {r) at r 



R, we have 
n\e\R 



J m (\U + e\R)N m+1 (\e\R) 

J m +i(\U + e\R)N m (\e\R)sign(e)sign(e + U ) 
ir\e\R 



J m (\U + e\R)J m+1 (\e\R) 



+J m +i(\U + e\R) J m (|e| J R)sign(e)sign(e + U 



(29) 



Using the asymtotics of the Bessel functions at large value of argument, we find the phase 
shift 5 m (e) = — arctan(z/ m / /x m ). Critical values of g are given by the solution of the equation 
Jm{gc) = at m > and J\ m \-i(g c ) = at m < 0. Taking into account the asymptotics 

21nx „ , , 2H(| m | _ 1)! 



N lx) 



\m\\ 



71 . . • - xH7T 

at x <C 1, we find for |e|i? 1, |e| C/q, and g = UqR close to g c 



*o(e) 



arctan 



fefl 



arctan 



7r( ej R) 2m+1 



2 2m m!(m - 1)! [(2m + 1) eR + 2m(^ - ^ c )] 



atra > 



(30) 



and 5_| m |(e) = 5| m |_i(e). If e < 0, which corresponds to electrons inside Fermi surface, and 
g < g c , then 5 m is always small. For g > g c , the phase shift 5 m (e) can be equal to ±7r/2 at 
some e < 0. That means the appearance at g = g c of the additional quasi-bound state on 
the Fermi surface. 

Calculation of the phase shift in the step-like potential was previously performed in 
Ref.— . However, the coefficients corresponding to fi m and u m , Eq. (l29|) . were found in Ref.— 
by matching the function u m (r) and its first derivative at r = R, instead of matching the 
functions u m (r) and d m (r). It is easy to check that the first derivative of u rn {r) is not a 
continuous function in the point r = R. As a consequence, the asymptotics of the induced 
charge density at large distances obtained in Ref.— is not correct. 



VI. AN INDUCED CHARGE 

Let us consider the induced charge Q>{r) outside of the radius r 3> R, 

POO 

Q>(r) = 2tt / dr'r'p ing {r') 

J r 

eNJ (g)J 1 (g)Rr 2 



7T 



de e' 



,Kl(er) + K (er)K 2 (er) - 2Kf 



er 



J 2 M + Jl{g){eR)nn 2 {eR) 



(31) 
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For |Jo(50l ^ (R/r)\n(r /R), we have 




(32) 



In the case \g — g c \ <C {R/r) ln(r/i?), we find with logarithmic accuracy 




(33) 



Since N = 4, then Q>(r)/e tends to the integer number iVsign(g — g c ) at g — > g c , having 
opposite sign for g < g c and g > g c . 

Let us discuss the induced charge Q<{r) inside of the radius r ^> R. Since the total 
induced charge Qtot = Q>{ r ) + Q<{ r ) is zero for the potential well at g less than the 
minimal g c , we have Q<{r) = —Q>{r) for such value of g. Note that Q to t is not zero 
for massless electron in graphene in the Coulomb potential Uc{r) = —Za/r even in the 
subcritical regime Za < 1/2, see Ref.— , due to zero mass of a particle and slow decreasing 
of a Coulomb potential at large distances. For g larger than the minimal g c , the total induced 
charge is already not equal to zero due to the effect similar to e + e~ pair production in the 
electric field of superheavy nucleus^i. In this case Q to t = eNM, where M is a number of g c 
less than g, so that Q < (r) = —Q>{r) + eNM. The quantity M is nothing but the number of 
the quasi-bound states at a given value of g, see discussion in SectiorJVl The explicit values 
of g c are given by zeros of the Bessel functions, as it is pointed out in SectiorJVl 

In Section IIVI and this section we have considered the contributions of the angular mo- 
menta m = and m = —1 in the Green's function (jSJ) to induced charge density and Q > (r) 
at large distances. Of course, the contributions of m > and m < — 1 are not zero though 
they are strongly suppressed by some power of R/r even in the vicinity of the corresponding 
critical points. However, M in Q tot = eNM includes numbers of g c coming from m > and 
m < — 1. 

In order to illustrate behavior of the induced charge in the vicinity of some critical point 
g c , it is convenient to represent Q>{r) and Q<(r) at r ^> R as follows: 



where F(g,r) is some continuous function of g. The dependence of this function on g at 
R/r = 0.1, obtained from Eq. (1311) in the vicinity of minimal value of g c , is shown in FigfS] 
(solid line), as well as its asymtotics, obtained with the use of Eq.(l32|) (dashed line). 



Q>(r) = eN[siga(g - g c ) + F(g,r)}, 



Q<(r) = eN[sign(g c -g)-F(g,r) + M], 



(34) 



10 



1 

.75 
^ 0.5 



-0.25 
-0 . 5 
-0 . 75 

2 2.2 2.4 2.6 2.8 
9 

FIG. 2: Dependence of the function F(g, r), denned in Eq. (|34|) . on g at R/r = 0.1 in the vicinity of 
a minimal value of g c . Exact result obtained from Eq. (|3ip is shown as a solid line, the asymtotics, 
obtained with the use of Eq. (|32p . as a dashed line. 

It is seen that the region, where Eq. (!32|) is not applicable, is very narrow. 

VII. CONCLUSION 

In this paper we have calculated the induced charge density generated by the potential 
well in graphene at large distances. Besides, we have obtained the induced charges outside 
of the radius r>i? and inside of this radius for subcritical and supercritical regimes. Small 
variation of the potential parameters drastically changes the induced charge density in the 
vicinity of the critical values of g. 

We are very grateful to O.P. Sushkov, G.G. Kirilin, and R.N. Lee for valuable discussions. 
The work was supported in part by RFBR grants 08-02-91969 and 09-02-00024. 
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